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Abstract. Porous-ceramic, thermal protection systems are used heavily in current reentry vehicles like 
the Orbiter, and they are currently being proposed for the next generation of US manned spacecraft, 
Orion. These systems insulate reentry critical components of a spacecraft against the intense thermal 
environments of atmospheric reentry. Additionally, these materials are highly exposed to space 
environment hazards like solid particle impacts. This paper discusses impact studies up to 10 km/s on 8 
lb/ft 3 alumina- fiber-enhanced- thermal-barrier (AETB8) tiles coated with a toughened-unipiece-fibrous- 
insulation/reaction-cured-glass layer (TUFI/RCG). A semi- empirical, first principals impact model that 
describes projectile dispersion is described that provides excellent agreement with observations over a 
broad range of impact velocities, obliquities and projectile materials. Model extensions to look at the 
implications of greater than 10 GPa equation of state is also discussed. Predicted penetration 
probabilities for a vehicle visiting the International Space Station is 60% lower for orbital debris and 
95% lower for meteoroids with this model compared to an energy scaled approach. 
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INTRO 

Porous-ceramic tiles insulate atmospheric 
reentry vehicles from reentry plasmas generated by 
atmospheric braking from orbital and exo-orbital 
velocities. Due to the necessity that these materials 
create a temperature gradient of approximately a 
thousand Kelvin over their thickness, it is 
important that these materials are near their pristine 
state prior to reentry. These tiles are also in general 
on exposed surfaces to space environment threats 
like orbital debris and meteoroids leaving a 
probability that these exposed surfaces will be 
below their prescribed values. Fortunately, owing 
to the typical small size of impact craters in these 
materials, the local flow fields over these craters 
afford some margin in thermal protection designs 
for these locally reduced performance values. 


The acceptability of a locally reduced 
thermal protection system is limited with the key 
limit being a direct impingement of the reentry 
plasma on spacecraft structure. For regions of the 
vehicle that are subjected to the most intense 
reentry environments, this limit of acceptability can 
also be realized even if the structure is not directly 
exposed as the structural elements can be heated 
above their safe operating condition. The 
remaining thickness of insulating material after an 
impact is then an important parameter describing 
the worthiness of the vehicle to reenter. As such, 
the depth of penetration is the principal observable 
required when testing the performance of these 
materials to the orbital debris and meteoroid 
environments. In the study reported here, NASA’s 
Hypervelocity Impact Test Facility (HITF) has 
impact tested these materials with projectiles 



typical of orbital debris and surrogate to 
meteoroids to determine the ballistic performance. 

EXPERIMENTAL PROCEDURE 

The tests performed are on target thermal 
protection materials typical of those used on the 
Orbiter and proposed for Orion. A density profile 
as a function of depth of a typical tile and 
substructure system is shown in Figure 2. From 
right to left, the figure shows the high density 
TUFI/RCG layer which tapers to the density of tile 
with its rear densification layer and the bonding 
pad and substrate to the extreme left. The 
TUFI/RCG layer provides a relatively high shock 
wave impedance material to push impacting 
materials to a higher pressure on impact. If the 
stresses are high enough to fragment/melt the 
projectile, the material decompresses extensively 
as it propagates in the highly porous tile resulting 
in scattered fragments and a diffusion of molten 
material. For threat particles too big to be arrested 
in the tile material, the densification and bonding 
layers’ increased density provide an increased 
arresting potential prior to structural panel impact. 

A pair of titanium powder enhanced 
orthogonal X-ray images of an off-normal impact 
crater is shown in Figure 1. The uppermost high- 
contrast plane in the views is the TUFI/RCG layer 
and the lowermost is the densification layer. The 
damage measurement normal to the TUFI/RCG 
layer is the maximum depth of penetration, and the 
damage measurements parallel to the TUFI/RCG 
layer are the widths. The product of the identified 
maximum depth of penetration and the tile density 
is the penetrated areal density value which can be 



Figure 2. Density profile experienced by projectile. 



Figure 1. Parallel (top) and perpendicular (bottom) 
cross section of an impact crater with cavity dimensions. 


used with the areal densities of the various layers 
of the thermal protection system. 

This effort has performed 50 shots where the 
main projectile cavity is within the tile and 
minimal bonding pad damage has occurred. 
Nylon™ (1.14 g/cm 3 ), aluminum 2017-T4 (2.796 
g/cm 3 ) and steel 440C (7.68 g/cm 3 ) projectiles have 
been considered. The penetrated cavity depth and 
widths are recorded for a variety of projectile sizes, 
impact obliquities, and impact speeds. The 
accuracy of the cavity measurements is limited to 
the resolution limit of the recording technique and 
to the limit of determining where cavities end. 
Using a fixed contrast ratio to define the cavity, the 
average measurement is used with the quadrature 
sum of the resolution limit and the variation about 
these averages for the measurement uncertainty. 
Along with the variable impact conditions, two 
areal densities of TUFI/RCG, 0.209±0.012 and 
0.158±0.009 g/cm 2 , have been considered. For 
partial penetration of the rear densification layer, 
the depth of the cavity includes an equivalent tile 
thickness of the penetrated densification layer. 

RESULTS AND DISCUSSION 

It has been observed that porous materials, after 
passage of a strong shock wave, jump to about the 
density of the solid material. This property of a 
shock wave manifests in a linear relationship 
between the shock wave velocity and the particle 
velocity without a constant offset or U = su where 
U, u and s are the shock wave velocity, the particle 
velocity and proportionality constant, respectively. 
Using this kinematic approximation in the shock 
wave frame, the Lagrangian momentum response 
to hydrodynamic and mechanical forces is 
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where x, M, p 0 , Y 0 and A are the shock wave 
penetrated depth, the propagating mass, the density 
and compression yield strength of pristine tile 
(-0.8 MPa) and the effective area of the bow shock 
wave, respectively. This equation is similar to that 
developed in Ref. 1; however, the propagating 
mass is considered here as the projectile, the 
TUFI/RCG layer and a conic of accreted tile 

M = ( m ~ + fnf - Po? ^ /a j n r p 2 + (2) 

and the effective area is also depth dependent. 

In Figure 3 the dependence of the ratio of half 
the cavity width to the penetrated depth for the 
normal impact records is shown as error ellipses. 
These expansion ratios are shown dependent on the 
ratio of an effective projectile areal density 
(m/ = 4/3 r p p p where r p and p p are the projectile 
radius and density, respectively) to TUFI/RCG 
areal density, m/. In this figure the cavity 
parameters are sorted by impact velocity bins with 
red, orange, green, blue and purple representing 4, 
7, 8, 9 and 9.5 km/s. The curves shown are 
empirical fits using a modified logistics function 
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where U t is the impact velocity and U m is a fit 
factor for this data set of 5.5 km/s. Nylon™ and 
steel projectiles have fit factors of 5.2 km/s and 2.6 
km/s, respectively. This function yields the 
required effective radius, r = co x + r p , which in 
turn yields the effective area. 

Transforming the penetrated depth in Eqn. 1 to 
the ratio of the effective to initial radius and 
normalizing by 1/2 p 0 U 0 2 results in 
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where f = r/r p , ( = U 2 /U 0 2 , ip = 2Y 0 /p 0 U 0 2 
and p = (m^ + rn^)/p 0 r p . The term U 0 is the 
initial shock wave velocity in the tile after release 
from the TUFI/RCG and derived from Eqn. 1 
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Figure 3. Dependence of aluminum impact radius to 
depth on projectile to TUFI/RCG for 4 (red), 7 (orange), 
8 (green), 9 (blue) and 9.5 (purple) km/s 
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where s T is the TUFI/RCG layer shock wave slope. 

The projectile and the swept tile mass are 
completely stopped when the dimensionless 
parameter for velocity goes to zero, thus the root of 
the function is the final non-dimensional radius. 
Transformation of these roots to projectile mass 
results in Figure 4. This figure represents the 
performance of this model relative to the collection 
of aluminum, Nylon™ and steel shots for all 
impact obliquities. To accommodate for the impact 
obliquity effects, a cosine of the impact obliquity 
term is used with the projectile density dependent 
powers of -2/3 for Nylon™, -8/3 for aluminum and 
-3 for steel. The dominate uncertainties in the 
abscissa are the penetrated areal density, rrif, and 
the expansion ratio which are related to the cavity. 

The model’s predicted ballistic performance for 
a 1.5” tile with a 0.209 g/cm 2 TUFI/RCG coating 
(heavyweight) is shown in Figure 5 in solid blue 
along with the ballistic performance for equal 
thickness tiles using 0.158 g/cm 2 (lightweight 
TUFI/RCG) and 0.044 g/cm 2 (RCG only) coatings 
in solid green and red, respectively. In all three 
cases, projectile sizes below the respective curves 
are not expected to fully penetrate the tile while 
projectile sizes above the curves are expected to 
penetrate the tile. Due to the longer duration at 
high stress states, the onset of fragmentation and 
melt with heavyweight TUFI/RCG is achieved at 
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Figure 4. Critical mass dependence on impact 
parameters with impact data represented by error 
ellipses and exact model agreement on the curve. 

lower impact velocities than with lightweight 
TUFI/RCG and RCG only layers. 

The shock wave compression is an important 
parameter to the ultimate ballistic performance of 
the tile. It has been shown that shock wave 
compression conditions transition between two 
different nearly constant values [2]. For lower 
shock wave strengths the material compresses to 
near the crystalline density, and at high shock wave 
strengths the material compresses to only a few 
times the initial density. This behavior has been 
incorporated in this impact model using a 



Figure 5. Ballistic limit curves for a 1.5” tile with 
heavyweight TUFI/RCG (blue) and equal thickness of 
lightweight TUFI/RCG (green) and RCG (red) surface 
coatings. Solid curve is compression to near-crystalline 
density and dashed is six-fold compression at 10 GPa. 


piecewise function for s on shock wave velocity 
where s is taken as 1.2 (six- fold compression) for 
shock wave velocities above 9 km/s and 1.055 
(~nineteen-fold compression) for velocities below, 
which is similar to the behavior shown in Ref. 4. 
The ballistic performance of the 1.5” tile with this 
behavior is shown as dashed curves in Figure 5 
with the same coloring. As can be seen, 
incorporating dissociation/ionization in the model 
predicts an increased performance. 

To determine the impact of these observations 
on spacecraft risk assessments, a comparison of the 
ballistic performance at a depth of 2.7 cm from this 
model and a traditional energy scaled model like 
that of Ref. 4 are used. Using these models with the 
solid particle environments for the International 
Space Station (ISS) indicates a 60% lower 
predicted probability of orbital debris penetrating 
to 2.7 cm in a tile with a lightweight TUFI/RCG 
coating. Due to the higher expected impact 
velocities of meteoroids, the predicted probability 
of penetration to the same depth is 95% lower. 

CONCLUSIONS 

Using the data obtained, a model using the 
material equation of state and strength properties of 
AETB8 tiles are used to develop a model that 
explains these findings and facilitates the further 
extrapolation to alternative tile configurations and 
impact conditions. Further work is necessary to 
develop analytical models of the projectile 
interaction with the discrete surface coatings of the 
tile and the obliquity effects, along with, further 
studies and work is necessary to extend the model 
to determine the impact of higher shock wave 
impedance at > 10 GPa on ballistic performance. 
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Referee Material 


Equation 1 represents a Langrangian form of the combination of the continuity, a simplified kinematic form 
of the equation of state and the momentum conservation equation. Starting with the momentum conservation 
equation in a frame at rest with the particle along the propagation direction 
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the left-hand side is the momentum changes experienced by the projectile and the target material, and the right- 
hand side is the sum of forces acting on the projectile. The mass within the control volume, M, includes the 
compressed material and the projectile. The kinematic form of the equation of state is U = su where s is the 
linear relationship between the particle velocity, u , and the shock wave velocity, U, assuming near zero sound 
speed. With this equation of state the hydrodynamic (differential pressure between shock wave front surface and 
vacuum surface) is given by: 


Fh — Po ^ 2 / 5 A. (A2) 

The mechanical (resistance to compression) force is given by: 

F m = Y 0 A. (A3) 

The terms p 0 , Y 0 and A are the initial density of the tile, the compressive yield strength of the tile and the 
projected area of the bow shock wave, respectively. The mass within the control volume is a summation of the 
mass of the projectile, TUFI/RCG layer and the accreted tile mass given by: 

M = ( m ~ + Wf- — ) n r p 2 + (A4) 

where m p and raf, are the effective areal density of the spherical projectile (4/3 p p r p where p p and r p are the 
projectile density and radius) and the TUFI areal density. The second term is the mass of the tile that has 
directly undergone accelerations by the expanding shock wave front from the projectile which is the product of 
the density of the tile and the volume of the cone created by the expansion where r p /co removes the volume that 
is extrapolated beyond the TUFI/RCG layer interface and w is from Eqn. 2. Transforming the left-hand side to 
Fagrangian coordinates and expanding yields: 
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where the control volume advances at the shock wave velocity. The first term in the expanded equation is the 
reaction of the control volume to the applied forces, and the second term is the momentum transfer to the 
i p 


accreted mass. The term : 
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transforms the accreted mass from the control volume to the frame at 


rest with the projectile accounting for the non-isentropic effects that occur within the shock front. Combining 
the expanded Fagrangian form of the momentum equation with the forces experienced by the advancing control 
volume on a per unit area basis yields 
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where the simplified tile equation of state property — = — is used. Using the projected radius equation 
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r = 60 x + r p to transform the independent variable to correspond to the initial condition and dividing through 
by 1/2 p 0 U 0 2 one arrives at 
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where U 0 is the velocity of the projectile as it breaks through the TUFI/RCG layer into the low density tile. 
Equation 4 is the non-dimensional form of A7 when r is transformed to = r/r p given by 
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where ^ = U 2 /U 0 2 ii = (m p + m T ) / (p 0 r p ) and p = Y 0 /(l/2 p 0 £/ 0 2 )- This first-order differential has the 
solution 
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The general solution for /is given by 
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where a t are the three roots of the equation 
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where k = (3 O) q) 1 ^ 3 . The combination of A9, A1 1 and A12 is the full analytic solution to Eqn. 4. For the case 
that k 3 is large, these roots are approximately given by 
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where i = 1, 2, 3. 

Performing an asymptotic expansion of the simplified form of All at the initial condition = 1 for k 
greater than one, the initial condition integral is given approximately by 
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Figure 6. Initial condition approximation (orange) relative to full numerical integral (black) as a function of 
the TUFI/RCG parameter and the error of estimate for the initial condition integral. 

The dependence of the 3 (s + 1) root of this approximate initial condition is shown in Figure 6 (orange curve) 
as a function of k relative to the full numerical solution of A1 1 (black curve). As can be seen in the Figure 6 the 
function is reasonably approximated to within a few percent to this simple linear relationship for values of k in 
the range of 2.5 to 4.5 which is typical of the aluminum impacts discussed in Figure 6 with the predominant 
value of A; being at about 3.5. Similarly in Figure 7, an asymptotic expansion at the condition <f » k > 1, yields 
the integral dependence 
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This approximation (orange curve) along with the full numerical solutions of All for A; (gray curves) equal to 
2.5, 3.5 and 4.5 (light, medium and dark, respectively) are shown along with the errors of estimate for k equal to 
2.5, 3.5 and 4.5 (light, medium and dark orange, respectively). Again it can be seen that this approximation is 
good for the typical > 10 expansions of Figure 7. Substituting these approximations into A9 one obtains 




Figure 7. Final condition approximation (orange) relative to full numerical integral for k equal to 2.5 (light 
gray), 3.5 (medium gray) and 4.5 (dark gray) as a function of the expansion and the error of estimate for k equal 
to 2.5 (light orange), 3.5 (medium orange) and 4.5 (dark orange) on the final condition integral. 
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Under the impact conditions considered here, i/j is much smaller than 1 allowing the numerator in the root to be 
approximately 1. Using the definition of f and solving for the projectile radius on the left-hand side yields 


o)x f Frr) rizL 1 1 ( s r o ycs+o 

(3 0 ) n) 1 / 3 ^ [9 V3 (3 0 ) n) 1 / 3 J Vs+1 1/2 p 0 t/ 0 2 / 


(A 18) 


where the dimensionless parameter definitions are substituted back with the exception of the dimensionless 
reduced mass. A Taylor expansion approximation is substituted for the small but not negligible constant. 
Converting to the mass of the impactor 
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Using Eqn. 5 for U 0 and substituting the definition of the reduced mass, //, for the effective projectile and 
TUFI/RCG areal densities the relationship between the initial condition velocity and the impact velocity Uj is 
given by 
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Substituting the approximation that the shock wave velocity is the particle velocity in the tile and the initial 
condition velocity into A19 yields 
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Noting that the value of s , is typically about 2 for a compression of an approximate 1.3 g/cm 3 TUFI/RCG layer 
one arrives at the simplified transcendental equation 
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which represents a 20% difference from the lead coefficient of 4/3 in Figure 4. 



